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Abstract 

Classically spinning B — 1 Skyrmions can be regarded as approximations to nu¬ 
cleons with quantised spin. Here, we investigate nucleon-nucleon scattering through 
numerical collisions of spinning Skyrmions. We identify the dineutron/diproton 
and dibaryon short-lived resonance states, and also the stable deuteron state. Our 
simulations lead to predictions for the polarisation states occurring in right angle 
scattering. 


1 Introduction 

The Skyrme model [T] is a nonlinear theory of pion helds that has topological soliton 
solutions. Encouragingly, Witten identified the Skyrme model as a low energy effective 
model of QCD [21E] . The conserved topological charge is interpreted as the baryon num¬ 
ber B, and the minimal energy static solutions for each integer B are called Skyrmions. 
They can be treated as rigid bodies, free to rotate in space and also in isospace (the 
three-dimensional space of the pion fields). When the rotational motion is quantised, 
the Skyrmions are models for nucleons and nuclei. In particular the 5 — 1 Skyrmions, 
quantised with spin and isospin half, model protons and neutrons. 
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A full quantum mechanical treatment of the interaction of two 5 = 1 Skyrmions is 
difficult. Each Skyrmion has three position coordinates and three orientational coordi¬ 
nates, so two-Skyrmion dynamics involves twelve coordinates jl]. A truncation to ten 
coordinates has been useful for modelling the deuteron bound state [5], but no proper 
discussion of two-nucleon scattering is possible with this truncation. 

An alternative is a classical approach, using the idea that a classically spinning 
Skyrmion is a reasonable model of a nucleon pm. The classical angular velocity of 
the Skyrmion is fixed to match the quantised spin and isospin of the nucleon. Scattering 
of non-spinning Skyrmions was first performed using an axially symmetric ansatz [8], 
the first full-field simulation was reported in P, and multi-charge Skyrmion scattering 
was considered in eg. The scattering of spinning Skyrmions representing nucleons was 
considered by Gisiger and Paranjape P, and they obtained analytical formulae for the 
scattering angle for large impact parameters. There has also been some recent work on 
multi-charge, purely isospinning Skyrmions na. 

A systematic, numerical investigation of Skyrmion scattering without spin, at mod¬ 
erate and small impact parameters m, confirmed that matter is exchanged between the 
Skyrmions when the Skyrmions are close, and also showed that a substantial rotation of 
the Skyrmion’s orientation can occur. Here, we numerically investigate the scattering of 
two classically spinning 5 = 1 Skyrmions to model proton-proton, neutron-neutron and 
proton-neutron scattering at various impact parameters, including head-on collisions. We 
are particularly interested in the change in the polarisation states of the nucleons when 
they scatter, and also in finding evidence, using our classical approximation, for the 
deuteron and for the known two-nucleon resonance states. 

The format of this paper is as follows. We first introduce the Skyrme model and its 
calibration. Then we discuss how to identify classically spinning Skyrmions as protons or 
neutrons. The next section presents the results of our numerical scattering of Skyrmions, 
and our identification of the dineutron/diproton resonances, the deuteron bound state, 
and the excited dibaryon state. The concluding section summarises and discusses our 
results. 


2 The Skyrme model 

The Skyrme model is defined by the Lagrangian density 

P2 1 2 p2 

£ = + ^TV((7 - h ), (1) 

where the Skyrme field U{t,x) is an SU(2)-valued scalar and 5^ = d/^UW is its su(2)- 
valued current. 5^., e and are parameters, which are fixed by comparison with ex¬ 
perimental data. Their values will be discussed later. It is convenient for us to work in 
dimensionless Skyrme units. One Skyrme length unit corresponds to ^ in inverse MeV, 
and one Skyrme energy unit corresponds to MeV. Conversion of inverse MeV to fm is 
as usual through h = 197.3 MeV fm. The dimensionless pion mass in Skyrme units is 

( 2 ) 
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In Skyrme units, the energy of a static field is 

E = j - ^Tv{[R„ Rj][Ri, R,]) + m^Tiih - U)^ d^x . 

It is often convenient, especially in numerical simulations, to express U in terms of a 
triplet of pion fields tt = (tti, 7 r 2 , tts) and an additional auxiliary field a as 

U{t, x) = a{t, x) I 2 + i7r(t, x) ■ t , (3) 

where r = (ti, r 2 ,T 3 ) are the three Pauli matrices, with the constraint cr^ + tt • tt = 1 so 
that U e SU( 2 ). 

At fixed time, U{t,x) is a map U : ^ SU( 2 ), where U I 2 sd spatial infinity. 

This boundary condition compactifies U { 00 } to S^. The group manifold of SU(2) is 
so a finite energy configuration U extends to a map U : ^ S^, and then belongs 

to a class of 7 r 3 ( 5 '^) = Z indexed by an integer S G Z, called the baryon number. B is 
also the degree of the map U which can be explicitly calculated as 

5 = y B{x)d^x = j £ijk^v{RiRjRk)d^x, (4) 

where B{x) is the baryon density. Skyrmions are the static field configurations of minimal 
energy for each value of B. In the figures we plot level-sets of baryon density B{x). 

It is well known that the B = 1 Skyrmion can be found using the so-called Hedgehog 
ansatz [T] 


Uh{x) = cos f{r) I 2 -\-f{r) X ■ T , (5) 

where r = |a;| and x is the radial unit vector, /(r) is a real radial profile function 
satisfying /(O) = 7r,/(oo) = 0. This produces a solution with rotationally symmetric 
energy and baryon density. Figure shows a H = 1 Hedgehog Skyrmion coloured as 
in [Zj. The centres of the white and black regions are where -|- 7r| = 0 and tts > 0, 
TTs < 0 respectively. The centres of the red, green and blue regions are where tts = 0 and 
tan“^ = 0, Y) ^ respectively. This is the colouring scheme used throughout this 
paper. 


3 O 

Side Top 

Figure 1\ B = 1 Hedgehog Skyrmion. 

An important feature of the Hedgehog ansatz is that an isorotation, U{x) AU (a?)Af, 
which acts as a rotation among the pion fields, is equivalent to a rotation in space. 
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U{x) U{D{A)x), where D{A)ij = |Tr(TjATjA'*') is the standard S0(3) rotation ma¬ 
trix corresponding to the SU(2) matrix A. This can be understood in terms of the 
coloured Skyrmions, where a spatial rotation can be ‘undone’ by a reordering of colours, 
i.e. an isorotation of the pion helds. A further use of the colouring is that it shows when 
two Hedgehog Skyrmions are in the attractive channel da. This is when the colours on 
the nearest, facing sides of two Skyrmions match. 

The B = 2 Skyrmion plays an important role in the dynamics of two 5 = 1 Skyrmions. 
This Skyrmion is toroidal, and is shown in figure 



Figure 2: B = 2 Skyrmion. 

To relate the Skyrme model to nuclear physics, Skyrmions should be quantised. As 
the Skyrmions are energy minima they can be treated as rigid bodies rotating in space 
and isospace. For general Skyrmions, semi-classical quantisation is the quantisation of 
the time-dependent space and isospace rotations A 2 and Ai in the rigid-body ansatz 

Uit, x) = Ar{t)Uo{D{A2{t))x)Ar{t )^, ( 6 ) 

where Uq{x) is a static solution. Ai{t) and A 2 (t) are SU( 2 ) matrices, and D{A2{t)) 
is, as before, the SO(3) matrix corresponding to A 2 {t). The spherical symmetry of the 
Hedgehog ansatz © implies that the rotational motion of the 5 = 1 Skyrmion only 
depends on the combined SU( 2 ) matrix A{t) = Ai{t)A2{t). 

Under rigid rotation the 5 = 1 Skyrmion has body-fixed angular momentum Lj = 
A(—Uj -|- bi) and body-fixed isospin angular momentum Ki = A(aj — fej), where aj = 
-iTr:{TjA\Ai) is the angular velocity in isospace, bj = i Ax{TjA 2 ^ 2 ) angular velocity 

in space and A is the moment of inertia M- From these one obtains the space-fixed 
angular momentum J* = —D{A 2 )JjLj and the “space-fixed” isospin angular momentum 
li = —D{Ai)ijKj. The space-fixed angular momenta correspond to physical spin and 
isospin, and these can be expressed purely in terms of A{t) and its time derivative. 
Nucleons have both spin and isospin half. A proton is an isospin-up state, and a neutron 
is an isospin-down state. 

Due to the spherical symmetry of the Hedgehog ansatz, Lj + = 0 {i = 1,2,3). 

Something similar is true for the 5 = 2 toroidal Skyrmion where L 3 -|- 2K^ = 0 |15j . 
We shall make use of these properties later. For an in-depth discussion of Skyrmion 
quantisation we point the reader to [isiiisiiniiii]. 
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2.1 Skyrmion calibration 

For the model to be relevant to nuclear physics it requires calibration. Adkins and 
Nappi [16] calibrated the quantised B = 1 Skyrmion against the masses of the spin | 
nucleons, the spin | delta resonances, and the pions, finding Ftt = 108 MeV, e = 4.84 
and m^r = 138 MeV (so m = 0.526). This calibration assumes that rigid rotation is a 
good approximation, but the surface fields of the delta are rotating close to the speed of 
light and pion radiation is very strong. Therefore, for the delta, the approximation is not 
reliable and this calibration is not valid for higher charge Skyrmions. Recently, Manton 
and Lau calibrated the Skyrme model against the states of Carbon-12 ra and found 
m = 0.7 to be optimal. For the B = 1 Skyrmion this leads to an energy (Skyrmion mass) 
M = 159.89, A = 55.89 and an rms matter radius = 1.006 in Skyrme units. The 

calibration requires = 117.5 MeV and e = 3.93. Planck’s constant in Skyrme units 
is always h = 2e^ so for this calibration h = 30.8. These are the values we use for our 
simulations and throughout this paper: 

m = 0.7 and h = 30.8. 

The resulting nucleon mass is the static B = 1 Skyrmion energy M plus the spin- 
energy contribution, M + which in physical units is 

With our calibration Mn = 1243 MeV, about 30% larger than the physical nucleon mass 
Mn = 939 MeV. Also the rms matter radius is = o.87 fm. 

The B = 2 toroidal Skyrmion can also be quantised as a rigid body. The two lowest- 
energy states are a spin 1, isospin 0 state representing the deuteron, and a spin 0, isospin 
1 state representing the low-energy resonance known as the diproton/dineutron. With 
our calibration, and following the quantisation procedure in ca. we find the mass for 
the deuteron to be 2324 MeV, which is similarly larger than the experimental value of 
1876 MeV. We also find the rms matter radius to be 1.08 fm, which is less than half 
the experimentally found rms radius of the deuteron ~ 2.14 fm, showing that rigid-body 
quantisation of the static 5 = 2 Skyrmion gives a too tightly bound representation of the 
loosely bound deuteron. A better representation of the deuteron in the Skyrme model 
was obtained in [5|. Using a Yang-Mills instanton ansatz for the Skyrme field [T9|, an 
additional, radial degree of freedom could be included, allowing the two 5 = 1 Skyrmions 
to separate. In the deuteron state, this radial degree of freedom oscillates with the zero 
point motion of a vibrational ground state. We will see later that modest oscillations of 
the radial degree of freedom are generated when spinning Skyrmions collide classically. 

Although this calibration is not optimal for the proton, neutron or deuteron, we have 
found that the processes discussed later are not sensitive to small changes of calibration. 


3 Spinning Skyrmions as protons or neutrons 

It was discovered by Adkins, Nappi and Witten [20| that the properly normalised wave- 
functions for the proton and neutron in spin-up and spin-down states, along the z-axis. 
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( 8 ) 

(9) 


=-{ai + ia2), 

TT 

= -{ao + ias), 
TT 


-(oo - ioa), 

TT 

-(oi - *02), 

TT 


where 2! = ao/2 + ^ai'rj is the SU(2) matrix that controls the 5 = 1 Skyrmion’s orientation 
and al + al + a 2 + al = 1. It has been proposed that these wavefunctions can be used 
to identify protons and neutrons as classically spinning Skyrmions PE]. This is because 
each wavefunction has maximal magnitude on a great circle in SU(2) and motion along 
the great circle is equivalent to a dyr rotation of the Skyrmion in space. The phase of the 
wavefunction changes by 27r around the circle, so the spin is 

The wavefunctions of p^ and are maximal when the Hedgehog Skyrmion (initially 
in its standard orientation) is simply rotated about the z-axis, whereas the wavefunctions 
of pf and n'*' are maximal when the Skyrmion is flipped over, and then rotated about the 
z-axis. 

A Hedgehog Skyrmion spinning with angular frequency uj about the z-axis is obtained 
by acting with the isorotation matrix A{t) = exp(^T3 ) = cos(y) + f sin(Y)T3- This 
isospinning Skyrmion is simultaneously spinning about its white-black axis in space. It 
has spin and isospin projections J3 = Xcu, I 3 = —Xcu, which are opposite. An anticlockwise 
spin, a; > 0, corresponds to an nX state and a clockwise spin, a; < 0, corresponds to a p'*' 
state. 

The isorotation matrix A(t){iT 2 ) flips the Skyrmion by tt about the p-axis, and rotates 
it with angular frequency u) about the z-axis. Such an isospinning Skyrmion has J^ = —Xuj 
and J3 = —Atu, so the spin and isospin projections are now equal. An anticlockwise spin, 
(u > 0, corresponds to an rA state and a clockwise spin, a; < 0, corresponds to a p^ state. 

The replacement A{t) —>■ A{t){iT 2 ) produces the spin flips rA —>■ u'*' and p'*' pf 
because D^z{A) —Dz‘i{A). This leads to the identification that a Skyrmion spinning 
anticlockwise about its white-black axis is always a neutron and a Skyrmion spinning 
clockwise about the same axis is always a proton, as shown in figure This is regardless 
of the orientation of the axis in space. 

To classically model a nucleon whose projected spin has magnitude | we require 
A|c<;| = so 



as h = 30.8 and A = 55.89. We use |a;| = 0.28 throughout this paper. With this angular 
velocity, the surface fields of the Skyrmion are not rotating close to the speed of light. 


l-t- I L I 't I I 


Figure 3: Identifying Skyrmions as neutrons and protons. 
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There are classical models of the two lowest-energy quantum states of the B = 2 
toroidal Skyrmion too. In the deuteron state, the torus spins around an axis orthogonal 
to the symmetry axis of the torus, with spin 1 . The colour orientation is arbitrary and 
unchanging in time, as the deuteron has isospin zero. In the dineutron/diproton state 
the torus is coloured as in figure and the colours circulate steadily, so that the isospin 
is 1. The black and white regions on the circular edge of the torus are fixed. The torus 
is not spatially rotating, as the dineutron/diproton has spin zero. The torus does not 
rotate about its symmetry axis in either state. If it did, the constraint L 3 + 2 K^ = 0 
would imply that the state had non-zero spin and non-zero isospin. 

4 Skyrmion scattering and its interpretation 

4.1 Spinning Skyrmion collisions 

We are interested in classically replicating nucleon scattering, namely proton-proton, 
neutron-neutron and proton-neutron scattering. To do this we identify protons and 
neutrons as spinning Skyrmions, as discussed above. Instead of colliding non-spinning 
Skyrmions as in some previous work m, we numerically collide spinning Skyrmions. We 
initiate a collision by forming two configurations Ui and U 2 out of spinning and translated 
Hedgehog Skyrmions, 

Uiit,x) = UH{Diit){x-X)), 

U2it,x) = UH{D2{t)ix + X)), (10) 

where Di{t) and T> 2 (t) are two SO(3) time-dependent rotation matrices and ±X are the 
locations of the Skyrmions. We choose X = (b,Y,0), where Y is large compared with 
the Skyrmion radius. Then, using the product ansatz, we create an initial configuration 
of two Skyrmions boosted towards each other, at velocities ±v, by 

U{t,x) = Ui{x,'y{y + vt),z)U 2 {x,'y{y -Yt),z) , (11) 

where 7 = 1 / y/1 — v^ is the usual Lorentz boost factor. 

Each Skyrmion is initially moving parallel to the y-axis. b is the impact parame¬ 
ter chosen to give the desired orbital angular momentum Ih = 2Mvb, where M is the 
Skyrmion mass. We define transverse polarisation to be where the Skyrmion’s white-black 
axis is aligned parallel to the z-axis and linear polarisation to be where the Skyrmion’s 
white-black axis is aligned parallel to the y-axis. 

We numerically evolve these initial configurations using the field equations obtained 
from the Skyrme Lagrangian density 0 . We use a finite difference leap-frog numerical 
algorithm on a suitably large lattice. Leap-frog was chosen because it is a symplectic 
integrator. 

All of the simulations are presented as videos at http://www.bristoltheory.org/ 
people/david.foster/skyrmevideos/; we urge the reader to view these videos. 

4.2 The Dineutron/Diproton 

There are no physical bound states in the 2 -proton or 2 -nentron systems. However, a 
2 -proton low energy resonance is observed in proton-proton collisions and is called the 
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diproton. It has spin J = 0 and isospin J = 1. The dineutron is the corresponding 2- 
neutron resonance. Neutron-neutron scattering is less easy to achieve than proton-proton 
scattering, but the dineutron resonance is observed as a decay product of the neutron- 
rich nucleus Beryllium-16, ^®Be where it is energetically favourable for ^®Be to decay 
to ^^Be by simultaneously ejecting two neutrons. The trajectories of the two detected 
neutrons trace back to the same point, where they form a state of lifetime ~ 10“^^ s. The 
^®Be and ^^Be ground states are both J = 0, hence the dineutron is a J = 0 state. The 
relative energy of the two neutrons is found to be ~ 0.25 MeV. So, assuming no other 
radiation, the neutrons separate with initial speed v = 0.02. 

There is no significant electrostatic repulsion between two neutrons, and since we 
have not included Coulomb effects in our Skyrme model, our simulations of Skyrmion 
scattering are likely to better model the dineutron rather than the diproton. So we 
collide two Skyrmions which are spinning as anti-aligned neutrons (illustrated in figure 
[^, giving a J = 0, / = 1,/s = —1 state, with initial speeds v = 0.02. The Skyrmions 
are either transversely polarised (parallel to the z-axis) or linearly polarised. Both types 


of scattering are 

shown in figure |4| 


Transverse 

polarisation 


• 

X 

^ 0 


t = to 

t — ti 

t = t2 

Linear 

polarisation 



X 

O 


t = to 

t — tl 

t = t2 

Figure 4: Level-sets of baryon density B{x) for the J = 

0 scattering of two neutrons. 


Three time-snapshots are shown here. The transverse polarisation scattering and the 
linear polarisation scattering videos are shown online. 

It is seen that no bound state is produced, but the B = 2 torus configuration is 
present for a short time. This torus forms in the (y, z)-plane and it does not spin, but it 
has isospin because the colouring continues to circulate at all times. It has the quantum 
numbers J = 0 and I = ff/s = —1 which is consistent with the short-lived dineutron 
state observed in the decay of ^®Be. For both polarisations, the black and white regions 
alternate around the circular edge of the torus. 


The scattering shown in figure |4] is for a head-on collision of two neutrons. For 
both polarisations, the outgoing Skyrmions are moving at right angles to the incom¬ 
ing Skyrmions, and they are spinning as neutrons too, because of isospin conservation. 
Similar right angle scattering was observed for non-spinning Skyrmions in m 

In the case of transverse polarisation the outgoing direction of motion is precisely 
the direction of the initial polarisation axis. There is another interesting phenomenon: 
the polarisation axis, the white-black axis, jumps by 90 degrees during the scattering, 
from being parallel to the z-axis to being parallel to the y-axis. This is possible because 
of the way the Skyrmions merge and break up, such that each outgoing Skyrmion is a 
combination of half of each incoming Skyrmion. 

From our numerical simulations we can observe the sense of the polarisation, before 
and after the collision. For example, the neutron incoming from the negative ^-direction 
(the left) has its polarisation in the positive z-direction, and the neutron that is outgoing 
in the positive z-direction has its polarisation in the negative y-direction. The result can 
be stated in a more invariant way, both for 2-neutron and 2-proton right angle scattering. 
Each particle, incoming or outgoing, has a momentum p and classical spin s. The two 
incoming particles and the two outgoing particles all have the same value for the vector 
p X s, so the scattering is accompanied by spin rotation. 

In the case of linear polarisation the outgoing direction of motion can be any direction 
at right angles to the y-axis, and is determined by the precise initial colour orientations. 
The scattering shown in figure has been selected to be from the ^-direction to the z- 
direction. Again there is a 90 degree jump in the polarisation axis, so that after scattering 
the neutrons are still linearly polarised. There is however a flip in the sense of the 
polarisation. In our example, for both incoming neutrons the polarisation is outwards, 
but for the outgoing neutrons it is inwards. Helicity, p • s, is therefore the same for 
incoming and outgoing particles, and this result is more general, holding in the diproton 
case too. 

Although the scattering here is classical, the results suggest that quantised 2-neutron 
and 2-proton scattering depends strongly on polarisations. For Skyrmions transversely 
polarised in the z-direction, the outgoing motion breaks the axial symmetry around the 
y-axis, because the outgoing particle motion is concentrated in the z-direction. Also, the 
outgoing polarisations are completely determined. It would be interesting to compare 
these classical simulations with physical proton-proton scattering in which the polarisa¬ 
tions of both incoming particles and both outgoing particles were measured. 

4.3 Deuteron formation 

The deuteron is a stable, two-baryon bound state comprised of a proton and a neutron. 
It has been understood as the quantised ground state of the B = 2 Skyrmion [laiain] 
with spin J = 1 and isospin / = 0. 

Here we are interested in classically modelling deuteron formation, so we numerically 
collide two spinning Skyrmions, representing a proton and neutron with transverse po¬ 
larisation parallel to the z-axis. The deuteron total angular momentum projected along 
the z-direction has three contributions, the individual spins Jp and of the proton and 
neutron, which are ±|, and the orbital angular momentum, an integer /, combining to 
give J3 = Jp + Jji + I [ 22 ] • Quantum models predict that the I = 0 state dominates, and 
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the deuteron spin is accounted for by the proton and neutron spins being aligned. We can 
replicate the deuteron state by colliding two Skyrmions head-on, where one Skyrmion is 
spinning as a neutron with Jn = \ and the other Skyrmion has the opposite orientation 
of the white-black axis and is spinning as a proton with dp = |. The impact parameter 
is zero, so / = 0. 

The deuteron has slightly lower energy than a proton plus a neutron, so we cannot 
use a simple energetic argument to choose an initial collision speed. Whatever the initial 
speed, some energy needs to be dissipated in order to truly produce a deuteron. It is 
therefore a considerable surprise that our collisions do produce a configuration similar to 
what is expected for a deuteron in a classical Skyrmion model. 

There is, however, a problem to be overcome. If we minimise the energy by choosing 
initial speed v = 0, and the initial relative orientations arbitrary, then averaged over time 
the Skyrmions repel. This is because the Skyrmions are spinning, and most of the time 
the colours do not match, producing a repulsion. Only briefly do the colours match, and 
the Skyrmions attract. Similar behaviour was observed for two spinning baby Skyrmions 
in 2-dimensions [23] • The repulsion of spinning Skyrmions initially at rest, which is not 
a Coulomb effect, is shown in figure 


X Q 9 9 3 ^ ^ 


t — to t — t — ^2 

Figure 5: The repulsion of a Skyrmion proton and a Skyrmion neutron with zero initial 
speed. The dynamics is shown online. 

Therefore, to model deuteron formation we need to collide the Skyrmions at a positive 
speed, and we choose v = 0.2. There is now enough kinetic energy to overcome the 
repulsive barrier. The dynamic 5 = 2 Skyrmion that forms is stable even with this extra 
energy. The collision along the y-axis of two Skyrmions spinning about the ; 2 -direction is 
shown in figure from two viewpoints. 
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Figure 6: Formation of a deuteron bound state in a spinning Skyrmion collision with 
J = 1, / = 0. This is shown dynamically online. 

As expected [ 21 ], the B = 2 torus forms in the (y,. 2 )-plane and spins about the 
axis, which is an axis orthogonal to the symmetry axis of the torus. Therefore there is no 
projection of spin or isospin along the symmetry axis, and the constraint L3 + 2K3 = 0 
is satished. 

The important feature of this collision is that the B = 1 Skyrmions do not scatter, but 
merge into a long-lived spinning 5 = 2 Skyrmion, with superimposed oscillations. The 
colours are not independently spinning. When the oscillation is at its maximal amplitude, 
two 5 = 1 Skyrmions can be identified. They represent two nucleons with parallel spins. 
This classical bound state models the physical, quantised deuteron na, with the correct 
spin, and zero isospin. The oscillation is similar to that of the 5 = 2 Skyrmion in 
its quantised ground state obtained using the instanton ansatz |3], but the amplitude of 
oscillation is somewhat greater, and the energy somewhat larger. The classical oscillation 
produces an oscillating matter radius, which is shown in figurej^as a function of time. The 
time-averaged matter radius is comparable with the experimental value for the deuteron. 

We have also attempted to simulate deuteron formation by colliding Skyrmions spin¬ 
ning as a neutron and proton, both with linear polarisation. This is shown in figure 
The process does not produce a classical bound state, even though the spin and isospin 
are correct for a deuteron. The initial state can be understood as Skyrmions approaching 
in the attractive channel with two white faces being closest. The Skyrmions scatter at 
right angles and the initial spin is converted to orbital angular momentum, the Skyrmions 
separating with some non-zero impact parameter. They also spin slowly about a coloured 
axis as they separate - this is inferred from a scattering of Skyrmions with larger initial 
spin. Skyrmions spinning like this should be interpreted as a superposition of a neutron 
and a proton. 
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t 

Figure 7: Skyrmion matter radius in the classical simulation of a neutron-proton collision, 
as a function of time t. The upper line is the experimental deuteron matter radius and 
the lower line is the matter radius of the static B = 2 toroidal Skyrmion. 


3 


t = tQ t — tl t = t2 

Figure 8: Head-on collision of linearly polarised proton and neutron with J = 1, I = 0. 
Note: at t = tl the torus is spinning about its black-black axis. The dynamics is shown 
online. 

Although the torus forms briefly in the linearly polarised case, it breaks apart. This is 
because Skyrmions in the attractive channel accelerate towards each other and the excess 
energy quickly breaks the torus up. In the transversely polarised case the excess energy 
is less and goes into the vibrational mode. This shows that polarisation has implications 
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for proton-neutron fusion. 


4.4 Excited Dibaryon 

Recently an excited proton-neutron state has been discovered [2S], called the dibaryon. 
It can be viewed as an excited spin 3 state of the spin 1 deuteron, with mass Moib ~ 2380 
MeV and quantum numbers J = 3, I = 0. 

To create the Skyrmion analogue requires colliding two Skyrmions each with kinetic 
energy T = |(MDib — 2 Mn); this is achieved with initial, relativistic speeds of v = 0.6. 
The J = 3, I = 0 state is replicated by colliding a Skyrmion spinning as a neutron and 
an oppositely orientated Skyrmion spinning as a proton, with the spins parallel and with 
orbital angular momentum 1 = 2, such that J 3 = ^ + ^+2. The orbital angular momentum 
is achieved with an impact parameter 6 = 0.32 at this speed. This impact parameter is 
the same order of magnitude as half the spatial width of the B = 2 toroidal Skyrmion. 
The scattering is shown in figure]^ with the initial motion parallel to the y-sods, and the 
proton and neutron polarisations parallel to the . 2 -axis. The dibaryon appears briefly in 
the form of the 5 = 2 toroidal Skyrmion spinning with J = 3 about the 2 ;-axis. The spin 
axis is at right angles to the torus’s symmetry axis, as for the deuteron. 
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Figure 9: Skyrmion scattering modelling proton-neutron scattering with J = 3 , / = 0 at 
the energy of the dibaryon resonance. The orbital angular momentum is / = 2, requiring 
an impact parameter b = 0.32. Note: at t = t 2 the Skyrmions are spinning about the 
red-cyan axis. The dynamical scattering is shown online. 


There is a small amount of excess energy and as a result the B = 2 Skyrmion almost 
breaks up into two rapidly spinning 5 = 1 Skyrmions moving back-to-back along the 
; 2 -axis. These have close to spin | each, but rapidly lose energy, presumably by pion 
radiation, although this is hard to see in the simulation. The spinning Skyrmions stop 
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and return to form a spinning B = 2 Skyrmion again. This replicates the experimental 
result that two pions are emitted, and a deuteron remains. 

The transient production of spin | particles matches other theoretical models of the 
dibaryon resonance involving two deltas with spin | [26]. Our simulation shows that 
(in the proton-neutron centre-of-mass frame) the deltas are emitted transversely to the 
line of collision, but parallel to the proton and neutron polarisation axes. Each delta is 
linearly polarised, and because its spin axis is an axis perpendicular to the white-black 
axis, it is in a superposition of states with different charges. This superposition is not 
surprising, because the average charge emerging along the positive or negative z-axis is 
half the proton charge. 

We have also attempted to produce a J = 3, / = 0 dibaryon state by colliding 
a proton Skyrmion and a neutron Skyrmion with the spin polarisations transverse but 
opposite, so that the spin angular momenta cancel. The orbital angular momentum is 
I = 3, which requires the same initial velocities, but an impact parameter b = 0.48. This 
scattering is shown in figure 



t = to t — t\ t — t2 

Figure 10: Modelling dibaryon production in proton-neutron scattering, with orbital 
angular momentum / = 3, and impact parameter b = 0.48. This is shown dynamically 
online. 

Here, the toroidal Skyrmion also appears but is not an excited deuteron, because it 
forms in the plane normal to the initial orbital angular momentum vector. In this case, 
one might anticipate the torus to be spinning about its symmetry axis, but it would then 
have non-zero isospin, which is not allowed. On closer inspection one observes that the 
B = 2 torus is not spinning, but each B = 1 Skyrmion can be interpreted as a wave 
propagating along the outer edge, such that the colours do not rotate. This curious 
dynamical configuration is not a dibaryon. 

In summary, dibaryon formation can be observed in classical Skyrmion scattering, but 
only if the initial polarisations of the proton and neutron are transverse and parallel. 
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5 Conclusion 


In this work we have classically scattered spinning Skyrmions to model nucleon-nucleon 
scattering. The initial spin angular momenta are fixed to be Interestingly, we have 
found short-lived states modelling the dineutron/diproton and the dibaryon. We have 
also found that in a collision with the quantum numbers of the deuteron, the Skyrmions 
do not scatter, but merge into a spinning and oscillating form of the B = 2 toroidal 
Skyrmion. So, even without quantisation, the Skyrme model usefully captures important 
physical features of two-nucleon dynamics. This is gratifying but a little surprising, 
especially when it is recalled that our long-lived deuteron state was made from two 
Skyrmions colliding with a larger kinetic energy than that needed to produce the short¬ 
lived dineutron/diproton. 

Our most valuable observations concern the polarisations. The classically spinning 
Skyrmions must be polarised along some axis, so we inevitably have complete polarisation 
information for the nucleons before and after scattering. Such complete information is 
hard to obtain experimentally, and also difficult to determine in more traditional nucleon- 
nucleon potential models. Classically we also fix the impact parameter. Several of our 
simulations lead to right angle scattering of the Skyrmions. Experimentally one cannot 
fix the impact parameter, but one can focus on scattering events where there is right 
angle scattering (in the centre-of-mass frame), and can hope that these are modelled by 
our collisions. For such events, our Skyrmion simulations make strong predictions for the 
polarisations. 

We have observed that in a 2-proton or 2-neutron head-on collision, with total angular 
momentum zero, if the particles are transversely polarised initially then they are trans¬ 
versely polarised finally. If they are linearly polarised initially they are linearly polarised 
finally. In both cases the polarisation axes rotate through a right angle. Moreover, there 
is a fixed plane containing the initial and final momenta and also the polarisations. The 
sense of the polarisations is determined by conservation of p x s for each particle in the 
transverse polarisation case, and conservation of p • s in the linear polarisation case. It 
would be very interesting if something similar was observed experimentally. 

In our simulations of a proton-neutron collision leading to the dibaryon, with total 
spin 3, we have seen the transient production of two linearly polarised deltas each with 
spin close to |. 

We intend to extend our numerical simulations to model nucleon-deuteron and deuteron- 
deuteron collisions, possibly forming Helium-3, Tritium or Hehum-4 bound states. The 
physical spins must be correctly modelled in the initial data, and the results will differ 
from those obtained with non-spinning Skyrmions. 
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